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We investigate the isospin chemical potential effect in the frame of SU{2) Nambu-Jona- 
Lasinio model. When the isospin chemical potential is less than the vacuum pion mass, 
the phase structure with two chiral phase transition lines does not happen due to C/a( 1) 
breaking of QCD. When the isospin chemical potential is larger than the vacuum pion 
mass, the ground state of the system is a Bose-Einstein condensate of charged pions. 
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1. Introduction 

It is generally believed that there exists a rich phase structure of Quantum Chro¬ 
modynamics (QCD) at finite temperature and baryon density, for instance, the 
deconfinement process from hadron gas to quark-gluon plasma, the transition from 
chiral symmetry breaking phase to the symmetry restoration phasID, and the color 
superconductivitj^at low temperature and high baryon density. Recently, the study 
on the QCD phase structure is extended to finite isospin density. The physical mo¬ 
tivation to study isospin spontaneous breaking and the corresponding pion super¬ 
fluidity is related to the investigation of compact stars, isospin asymmetric nuclear 
matter and heavy ion collisions at intermediate energies. While there is not yet pre¬ 
cise lattice result at finite baryon density due to the Fermion sign probleni^, it is in 
principle no problem to do lattice simulation at finite isospin densitj®. It is found^ 
that the critical isospin chemical potential for pion condensation is about the pion 
mass in the vacuum, /rj ~ TOtt. The QCD phase structure at finite isospin density is 
also investigated in many low energy effective models, such as chiral perturbation 
theorcfifladder OClS. random matrix methocP. stron. couolin. lattice Ocflini 
and Nambu-Jona-Lasinio (NJL) modelllllliail*!. 

One of the models that enables us to see directly how the dynamic mechanisms 

rm 

of chiral sym metry breaking and restoration operate is the NJL models applied 
to quarki^. Within this model, one can obtai n the hadronic mass spectrum and 
the stat ic proper ties of mesons remarkably weli ^^l^^ , and the chiral phase transi¬ 
tion lin J^ plane is very close to the one calculated with lattice 
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qclISOI. 

Recently, this mod el is used t o investigate the color superconductivity 
at moderate baryon densit}i^ ^ l ‘^^ l ^^ l ^^ l ^^ l In this letter we report our results on 
the phase structure of the general SU{2) NJL model with 17 a (1) breaking term 
at finite isospin chemical potential. At small isospin chemical pote ntia l, we argue 
that there is only one chiral phase transition line in the T — /is planJS. At isospin 
chemical potential larger than vacuum pion ma ss, the ground state of the system is 
a Bose-Einstein condensate of charged pionJl^. 


2. Mean Field Approximation at Finite Isospin Chemical Potential 

We start with the flavor SU{2) NJL model defined by 

£ = - mo + M7 o)'0 + , (1) 

where mo is the current quark mass, fj, the chemical potential matrix in flavor space, 

^ = diaginu, /id) = diag ^ 

with g,B and /r/ being the baryon and isospin chemical potential, respectively, and 
the interaction part includes^ the normal four Fermion couplings corresponding to 
scalar mesons cr, ao,a+ and a_ and pseudoscalar mesons r;/, 7ro,7r+ and 7r_ excita¬ 
tions, and the’t- HoofipSI 

determinant term for Ua{^) breaking, 

G ■ ^. - - K - 

Gint = -^ X! + (V'* 75 'ra' 0 )^] + Y [det ?/;(l -|- 75 )/)’ -b det ?/;(l - 75 ) 1 /'] 

a— 0 

= ^{G + K) +^{G- K) • (3) 

At zero isospin chemical potential jii = 0, for AT = 0 and mo = 0 the Lagrangian is 
invariant under C/b(1) 17a(1) 0 SUv{2) 0 SUa{‘2) transformations, but for K ^ 

0, the symmetry is reduced to U b{^) ^ SUv{2) ^ SUa{2.) and the C/a(1) breaking 
leads to a and a mass splitting and tt and 7 / mass splitting. If G = AT, we come 
back to the standard NJL mode^^l with only a, ttq, 7 r+ and 7 r_ mesons. 

We introduce the quark condensates 

Gu = {uu) , Gd = (dd) , (4) 


and the pion condensate 


TT - - 1 - 


(5) 


where we have chosen the pion condensate to be real. The quark condensate and 
pion condensate are, respectively, the order parameters of chiral phase transition 
and pion superfluidity. In mean field approximation, the thermodynamic potential 
of the system is 


n — G((T^ + (T^) + 2KGuGd + 




V 


2 


( 6 ) 
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where is the inverse of the mean field quark propagator, in momentum space 
it reads 


f + flu'lo - rriu ijziG + K)Tr \ 
\ ij5{G + K)TT fldlo - md) 


with the effective quark masses 


(7) 


TO„ = mo — 2Gct„ — 2Kud , md = mo — 2Gad — 2Kuu ■ ( 8 ) 


The condensates (t„ , Ud and tt as functions of temperature and baryon and isospin 
chemical potentials are determined by the minimum thermodynamic potential, 


an _ _ n _ n 

Tv Tv Tv ^ • 

UGu 


(9) 


It is easy to see from the chemical potential matrix and the quark propagator matrix 
that for /rs = 0 or /r/ = 0 the gap equations for (t„ and Ud are symmetric, and one 
has 


= CTd = a 12 , m„ = TOd = TO = Too - (G + K)a . (10) 


3. Chiral Phase Structure at fip < 


We now consider the QCD phase structure below the minimum isospin chemical 
potential /if = m^ for pion superfluidity. Since the pion condensate is zero, there is 
only chiral phase structure in this region. 

It is well known that in absence of isospin chemical potential the chiral conden¬ 
sate jumps down suddenly at a critical baryon chemical potential, which indicates 
a first order phase transition. At both finite baryon and isospin chemical potential, 
in principle we have cr„ ^ ad, they may jump down at the same critical baryon 
chemical potential or at two different critical points. In the case K = 0, the two gap 
equations for (t„ and ay decouple and the two critical points do exist, and therefore, 
there are two chiral phase transition lines in T — ft b plane at fixed isospin chemical 
potential, the interval between the two points is just A/ig = 3/i/. What is the effect 
of the Ua{1) breaking term on the QCD phase structure? If A" 7 ^ 0, the two gap 
equations couple to each other and tend to a single phase transition. Especially in 
chiral limit, one can clearly see that when one of the quark condensates becomes 
zero, the other one is forced to be zero for any coupling constants G and AT 7 ^ 0. 
Therefore, there is only one chiral phase transition line at any A" 7 ^ 0 in chiral 
limit. In real world with small current quark mass, numerical solution of the gap 
equations shows that if there exists a structure of two chiral phase transition lines 
depends the quantity 


G^Y 

G + K 2 V G + Kj 


( 11 ) 


In real world there are four parameters in the NJL model, the current quark 
mass Too, the three-momentum cutoff A, and the two coupling constants G and 
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K. Among them mo, A and the combination G K can be determined by fitting 
the chiral condensate a, the pion mass and the pion decay constant /,r in the 
vacuum. For a = 2(—241.5 MeV)^, = 140.2 MeV, and /,r = 92.6 MeV, one has 
mo = 6 MeV, A = 590 MeV, and (G + K)t\?/2 = 2.435. To determine the two 
coupling constants separately or the ratio a, one needs to know the r]i- or a-meson 
properties in the vacuum. In the self-consistent mean field appr oxim ation, the r]' 
mass mrii can be calculated in random phase approximation(RPA^|^ 

l_(G-iF)nps(fco = m^-,k = O) = 0 (12) 


with polarization function IIps defined as 
nps(/co,k) = -i 

The quantity a can be written as 


d'^p 


Tr 


lbSm.f{p + k)^bSmf{jp) 


(13) 




G 


Ki2f 


d^k 


Ek 


(27r)3 El-m'^,14 


(14) 


where Ek = Vk^ -|- m^. Combining with the above known parameters mo, A and 
G -I- AT obtained and choosing m^, = 958 MeV, we have a = 0.29 which is much 
larger than the critical value 0.11 for the two-line structure at /i/ = 60 Me’V^I. in 
fact, for a wide mass region 540 MeV < rurji < 1190 MeV, we have a > 0.11, and 
there is no two-line structure. 

To fully answer the question if the two-line structure exists before the pion 
condensation happens, we consider the limit pi = m^. In real world with mo yf 0 and 
AT 0, with increasing coupling constant K or the ratio a the two lines approach 
to each other and finally coincide at about a = 0.21 which is still less than the 
value 0.29 calculated by fitting m,,/ = 958 MeV. Therefore, the two-line structure 
disappears if we choose m^/ = 958 MeV. In fact, for 720 MeV < m,,/ < 1140 MeV, 
we have a > 0.21 and the two phase transition lines are cancelled in this wide 
mass region. When the ry/-meson mass is outside this region, the C7 a( 1) breaking 
term is not strong enough to cancel the two-line structure, but the two lines are 
already very close to each other. Considering the relation between the r]' mass and 
the number of flavors, m^, oc Vy, one can estimate the rj' mass in the case of two 
flavors m,,' « 780MeV which is in the region where a > 0.21. 


4. Pion Superfluidity at p/ > 

In this section we concentrate on the case /ip =0. The gap equations for chiral 
condensate a = au + Ud and pion condensate tt are derived as following 

^ -^(2/(A+)-1) 


a = 6m 




I -h 12iJ7r 


d^k 

(27r)3 


^(2/(A-)-l) + ^(2/(A+)-l) 


= 0 . (15) 
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where H = {G + K)/2,Ef = Ek ± fii/2,E^ = ^. The isospin 
density nj can be derived from thermodynamic potential 


nj = 3 


d^k 


(2/(E+) - 1) - (2/(E-) - 1) 

-t^TT dlz-jY 


(16) 


Numerical solution of the gap equations at zero temperature is shown in Fig. ^ 
We found that when /j,/ < m,r, the ground state is the same as the vacuum and 
the isospin density is zero, while pion condensate becomes nonzero when /i/ > 
and isospin density becomes nonzero. The phase transition is of second order. The 
critical isospin chemical potential /ij is exactly the vacuum pion mass which 
can be proved analytically. The critical isospin chemical potential is determined by 

1 - 12iL f - = 0 . (17) 


(27r)3 Ek El - {ijiff /4 


which is just the same as the mass equation for pion in vacuum. 

The collective meson excitations can be calculated in the framework of random 
phase approximation. Different from normal RPA approach in vacuum and at finite 
temperature, the meson modes in the pion superffuid will mix with each other and 
the dispersion relations of each eigen modes are determined by 


det(l- 2£rn(fco,k)) = 0 


(18) 


where the polarization function matrix n(fco,k) is defined as 


1 - 2iJn(fco,k) 


/1 - 2i7n,, 

-2i7ncr7r+ 

-277n+^_ 

-277n++„ 


1 - 2770,7^,7+ 

-277n,r_^,r_ 

-277n+++o 

-2Hn^_, 

-277n7r_,r^ 

1 - 277n+_+_ 

-277n,r_+n 

\ -2i7n7r(,cr 

-2770,7077+ 

-277n+o+_ 

1 - 277n+„+„ 


The polarization functions are defined as 

nMM'(^0ik)=i , (20) 

with the vertexes Tm and F^ 


Fm 


1 M = a ( 1 M = a 

fT+75 M = TT+ p* _ J *'^-75 M = TT+ 

M = 7r_ ’ ^ I iT +75 M = 7r_ 

iT375 M = TTo [ *T375 M = ttq 


( 21 ) 


where t± = (n ± iT2)lV^. From detailed calculation one has 

^cr7r± — llTr^CT OC (TTT, Il7,.^7,._ — Il7r_7r+ OC TT , — 11/7,.^ — 0(7 — ^i^it)- (22) 

Thus ttq does not mix with other modes, while ct, 7r_|_,7r_ mix with each other. 
However, mixing between a and tt± is only strong near the phase transition point 
and can be neglect at large /r/. The meson mass spectrum as a function of fij is 
shown in Fig. |21 We found that the full calculation (solid line) is consistent with 
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Fig. 1. The order parameters cr and tt, scaled by the vacuum value (tq, and isospin density, scaled 
by the value of baryon density in normal nuclear matter. 

chiral perturbation theory, and the mixing between cr and 7r± plays a very important 
role. On the other hand, in the superfluid phase with tt 0, det(l — 2iLn(0,0)) = 0 
is always satisfied, which indicates that Goldstone’s theorem is satisfied in the whole 
superfluid phase at RPA level. 



Fig. 2. Meson masses as a function of isospin chemical potential. The solid line is full calculation 
and the dashed line is obtained by neglecting a — n± mixing. 


5. Summary 

We have investigated the two flavor NJL model with 17a(1) breaking term at hnite 
isospin chemical potential. The main conclusions are: 

1) The two chiral phase transition lines in the T — plane predicated by the NJL 
model without Ua{^) breaking term are cancelled by the strong 17 a( 1) breaking 
term at low isospin chemical potential. Therefore, in relativistic heavy ion collisions 
where the typical value is much less than m.^., it looks impossible to realize the 
two-line structure. 

2) The critical isospin chemical potential for pion condensation in NJL model is 
exactly the pion mass in the vacuum, /ij = independent of the model param¬ 
eters, the regularization scheme, the 17a(1) breaking term. When fij exceeds m^, 
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The ground state of the system is in a Bose-Einstein condensate of charged pions. 
In this pion superfluid, the meson modes cr, 7r+, 7r_ mixed with each other and there 
is a Goldstone mode due to the spontaneous breaking of Ui{l) symmetry. 
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